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The importance of thermoelastic damping as a fundamental dissipation mechanism for small- 
scale mechanical resonators is evaluated in light of recent efforts to design high-Q micrometer- and 
nanometer-scale electro-mechanical systems (MEMS and NEMS). The equations of linear thermoe- 
lasticity are used to give a simple derivation for thermoelastic damping of small flexural vibrations 
in thin beams. It is shown that Zener's well-known approximation by a Lorentzian with a single 
thermal relaxation time slightly deviates from the exact expression. 
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I. INTRODUCTION 

Micro-electro-mechanical systems (MEMS) and more 
recently nano-electro-mechanical systems (NEMS) are 
being developed aggressively for a variety of applications 
as well as for accessing new regimes of basic experimental 
research. Among the different applications envisioned for 
MEMS and NEMS are ultrafast and high precision actu- 
ators, sensors (such as accelerometers, bolometers, mag- 
netometers, and calorimeters), and narrowband high- 
frequency mechanical filters, all with compact and low- 
power designs that can be fully integrated with modern 
semiconductor electronics. Experimentally, it is hoped 
that NEMS will open the door to the investigation of 
new regimes of phonon-mediated processes as well as ' 
quantum behavior of mesoscopic mechanical systems.l 

For all these pursuits it is desired to design and con- 
struct systems with very little loss of energy or very 
high quality factors Q. Unfortunately, it has been con- 
sistently observed that the quality factors of resonators 
decrease with size significantly — even when made from 
pure single-crystal materials. It is therefore of great im- 
portance to understand the dominating energy dissipa- 
tion mechanisms in mechanical resonators when one ap- 
proaches submicron scales. Within these mechanisms one 
would like to identify those which are fundamental and 
always impose an upper limit on the quality factor and 
those which might be eliminated through improved de- 
sign and fabrication. 

In this work we examine the importance of the pro- 
cess of thermoelastic damping as a dissipation mecha- 
nism in MEMS and NEMS. Although there has been 
some cUsQussion of different dissipation mechanisms in 
MEMSuTa very few authors have addressed the question 
of thermoelastic damping. Roszharta observed thermoe- 
lastic damping in single-crystal silicon miojo-resonators 
at room temperature; and Yasumura et aln recently re- 
ported thermoelastic damping in micro-resonators an or- 
der of magnitude smaller than Roszhart's in silicon ni- 
tride, also at room temperature. The question arises 
whether one should be surprised by such observations or 



whether one should expect to see thermoelastic damping 
at these length scales. We establish in this paper that 
as long as the system remains in the regime of "diffusive 
thermal phonons" (to be discussed below) the latter is 
the case for two basic reasons: 

1. The strength of the damping caused by thermoelas- 
tic coupling is independent of geometry. It depends 
only on the thermodynamic properties of the ma- 
terial as a function of temperature. 

2. In the case of flexural vibrations of thin beams the 
position of peak damping as a function of frequency 
depends on the dimensions of the beam. There- 
fore, even though the normal frequencies of the res- 
onators increase as they become smaller so does the 
frequency at which peak damping occurrs. 

These effects conspire together to maintain the rele- 
vance of thermoelastic damping all the way down to the 
nanometer scale. 

In the next section we describe the process of thermoe- 
lastic damping, review some of the relevant literature, 
and present the outline of this paper. 



II. THE PROCESS OF THERMOELASTIC 
DAMPING 

Acoustic modes — such as a sound wave traveling 
through an infinitely large elastic material or a nor- 
mal mode of vibration of an elastic resonator of finite 
geometry — will experience damping due to their nonlin- 
ear interaction with a surrounding bath of thermally- 
excited elastic modes, or phonons. If the mean free path 
of these thermal phonons is much smaller than the wave- 
length of the acoustic mode then sufficient thermaliza- 
tion occurrs on the scale of interest. It is then possible 
to define a temperature locally, even when the system 
is not in a state of thermal equilibrium. Equivalently, if 
the relaxation rate of the phonon distribution to a local 
Bose-Einstein distribution is much faster than the fre- 
quency of the acoustic mode then one has a well-defined 
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temperature field, and there is no need to treat the ther- 
mal phonons as individual excitations. In this regime, 
in which thermal phonons are said to be "diffusive," the 
complicated interaction between the acoustic mode and 
the thermal phonon bath in an isotropic solid is captured 
by a single macroscopic parameter — the material's ther- 
mal expansion coefficient 
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which couples changes of length to changes of tempera- 
ture. Note that here and throughout the paper we make 
use of the linear coefficient of thermal expansion whose 
value is one third of the volumetric coefficient of thermal 
expansion. 

When an elastic solid is set in motion it is taken out 
of equilibrium, having an excess of kinetic and poten- 
tial energy. In an isothermal and perfectly-linear elastic 
solid such a non-equilibrium state can exist forever. In 
a thermoelastic solid the coupling of the strain field to a 
temperature field provides an energy dissipation mecha- 
nism which allows the system to relax back to equilib- 
rium. Relaxation of the thermoelastic solid is achieved 
through the irreversible flow of heat driven by local tem- 
perature gradients that through the coupling accompany 
the strain field. This process of energy dissipation, called 
thermoelastic damping, is a fundamental one. As long 
as a — which acts as a coupling constant — is non-zero, 
thermoelastic damping introduces an upper limit to the 
quality factor of even the most perfectly designed and 
constructed resonator. 

The first to realize that thermoelastic relaxation may 
be a significant source of damping in mechanical res- 
onators was Zener, who in a series of papersEJ in the 
1930's developed a general theory of thermoelastic damp- 
ing. Zener treated the problem in the frameworl^f his 
so-called "standard model" of the anelastic solidlJ and 
showed that the damping behavior can be approximated 
very well by a single relaxation peak with a characteris- 
tic relaxation time. This relaxation time corresponds to 
the thermal diffusion time across the width h of the beam, 
which is proportional to where x is the solid's ther- 

mal diffusivity [to be defined in Sec. |v| Eq. (|l3|)]. We 
shall review Zener's theory in section pll[ 

It was only two decades later that other researchers be- 
gan reexamining the question of thermoelastic damping 
by seeking exact solutions to the coupled equations of lin- 
ear thermoelasticity in various geometries. The solutions 
for ppoBagating plane waves in an infinite thermoelastic 
solidli^ t3 showed that the two transverse modes, which 
propagate through the solid without causing any local 
volume changes, do not couple to the temperature field 
and hence suffer no damping. The longitudinal mode, on 
the other hand, does couple to the temperature field and 
an exact expression has been obtained for its attenua- 
tion and dispersion. Solutions for thermoelastic Rayleigh 
waves (two-dimensional surface waves on a semi-infinite 



solid) immediately foUowedjli but progress on thermoe- 
lastic solids with finite geometries came much later due to 
the well-known difficulty of solving even the equatipns of 
linear elasticity with finite boundary conditionsO Nev- 
ertheless, analytic solutions now exist, at least in terms 
of the dispersion relations, for thermoelastic waves in an 
infinite thin plate,E3 and longitudiDaljjjjaves in infinite 
rods with circular cross sections. tSEarJ Other geome- 
tries, such as beams of rectangular cross sections, have 
been too difficult to solve analytically. To treat such 
problems one generally needs to use approximate theo- 
ries. 

MEMS and NEMS resonators generally contain ele- 
ments which vibrate in either torsional or flexural modes. 
Because of the way one fabricates such devices, rectangu- 
lar cross sections often turn out to be the most relevant. 
Pure torsional modes of rectangular beams involve no lo- 
cal volume changes and therefore, just as for transverse 
waves in the bulk, they do not suffer any thermoelastic 
losses. For fiexural vibrations of thin rectangular beams 
one may use Zener'alJ approximate expression [Sec. Ill 
Eq. (|^)] for thermoelastic damping. This appears to be 
the general practice.™ On the other hand, one may try to 
seek exact solutions of the thermoelastic equations for the 
case of a thin beam under flexure. Landau and LifshitzcJ 
have provided an exact expression for the attenuation co- 
efficient of such vibrations without, however, giving a rig^ 
orous derivation and solution of the equations. ShichEj 
investigated the thermoelastic beam equations in a more 
general context while performing a dynamic instability 
analysis of vertically-standing cantilevers. He solved the 
equations, giving a plot of thermoelastic damping 
in thin rectangular beams, but fell short of providing an 
analytic expression for . 



To remedy this state of affairs we give, in section IV, a 



simple derivation of the approximate thermoelastic equa- 
tions for a thin beam under flexure, and then solve these 
equations in section ^ to arrive at an exact expression 
[Sec. 0Eq. (29)] for thermoelastic damping in thin rect- 
angular beams. Although Zener's approximation is good 
for many situations we suggest that in the future the ex- 
act expression given here be used instead. In section VI 
we discuss the experimental implication of our results 
showing the characteristic damping curves expected for 
GaAs and Silicon, which are typical materials used in 
the fabrication of micrometer- and nanometer-scale res- 
onators. In section VII we say a few words about ther- 



moelastic damping of l ongitu dinal waves in MEMS and 



NEMS, and in section VIII we conclude by discussing 



the validity of the theory which we present here. In the 
appendix we give a careful analysis comparing the ex- 
act expression for thermoelastic damping in thin beams 
with the approximate Lorentzian behavior, predicted by 
Zener. 
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III. ZENER'S STANDARD MODEL OF THE 
LINEAR ANELASTIC SOLID 

Zener's standard model of anelasticitj0 is based on 
an extension of Hooke's law to the most general linear 
homogeneous equation involving stress tr, strain e, and 
their first time derivatives: 



(2) 



The physical interpretation of the three parameters of the 
model is straightforward. When the strain is kept con- 
stant the stress relaxes exponentially with a relaxation 
time r£. Similarly, To- is the strain relaxation time when 
the stress is kept constant. is the value of the per- 
tinent elastic modulus after all relaxation has occurred. 
Mu — Mii{Ta /Tf) is the unrclaxed value of the elastic 
modulus. 

Under periodic dynamical conditions 



a{t) = (Toe'' 



(3) 



the stress and strain amplitudes are related by a 
frequency-dependent complex elastic modulus. The dis- 
sipation, or "internal friction," in the solid is defined 
as the fraction of energy lost per radian of vibration. If 
it is small, the dissipation is equal to the ratio of the 
imaginary and real parts of the complex modulus, giving 



where r — ^/rvTe, and 
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is a dimensionless quantity called the "relaxation 
strength" of the modulus. 

Thus, the dissipation exhibits a Lorentzian behavior as 
a function of lot with a maximum value of AAf/2 when 
ujT = 1. Dissipation peaks of this form, generally called 
"Debye peaks," are quite ubiquitous as to be expected 
from their prediction by such a nai've model. They occur 
for many different relaxation mechanisms such as point- 
defect relaxation ("Snoek peaks"), defect pair reorien- 
tation ("Zener peaks"), dislocation relaxation ("Bordoni 
peaks"), grain |boundary relaxation, and of course ther- 
mal relaxation.c2l In many of these examples there is not 
just a single relaxation time r and therefore one sees mul- 
tiple or broadened Debye peaks. 

One can understand qualitatively why there is a peak 
in dissipation when lot = 1 in the following way: If the 
frequency of vibration to is much smaller than the ef- 
fective relaxation rate 1/r of the solid then the system 
remains essentially in equilibrium and very little energy 
is dissipated. If the vibration frequency is much larger 
than the effective relaxation rate lo ^ l/r, the system 
has no time to relax and again very little energy is dissi- 
pated. It is only when the vibration frequency is on the 



order of the system's effective relaxation rate that appre- 
ciable dissipation occurs. The full picture may be more 
complicated, however, because in some cases r itself can 
depend on lo. 

In the case of a thermoelastic solid the relaxation 
strength (^) to be considered is that of Young's modulus 



Egd — E 

E 



(6) 



whose value is known from basic thermodynamics. Here 
Ead is the unrclaxed, or adiabatic, value of Young's mod- 
ulus and E is its relaxed, or isothermal, value. Cp is the 
heat capacity per unit volume at constant pressure, or 
stress, but replacing it by the heat capacity at con- 
stant volume, or strain, will introduce an error in A^; 
which is only on the order of Ae^- Since , and 
therefore also A^;, are assumed small such an error is 
negligible. For similar considerations, no harm is done 
by replaciiag \/EadE in the denominator of (||) by E. 

ZeneiUj'EJ calculated the thermal relaxation times as- 
sociated with different transverse thermal modes for a 
thin beam under flexure. He showed that for rectangu- 
lar beams approximately 98.6% of the relaxation occurs 
through the first mode whose relaxation time is 



(7) 



where x is the thermal diffusivity of the solid and b is the 
width of the beam. Only a very small error is therefore 
made by consideringthe vibrating thin beam as having 
a simple relaxation (j^ with a single relaxation time r^. 
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Nevertheless, we show below that for the simple geom- 
etry of a rectangular beam — which is the most relevant 
for many current MEMS and NEMS designs — such an 
expansion in transverse thermal eigenmodes is unneces- 
sary and an exact expression for thermoelastic damping 
can easily be obtained. 



IV. THERMOELASTIC EQUATIONS OF A THIN 
BEAM 

We consider small flexural displacements of a thin elas- 
tic beam of length L and rectangular cross section of di- 
mensions b X c. We define the x-axis along the axis of 
the beam and the y- and z-axes parallel to the surfaces 
of dimensions b and c, respectively. In equilibrium, the 
beam is unstrained, unstressed, and at temperature Tq 
everywhere. Departure of the beam from equilibrium is 
described by a displacement field Ui (i = x,y,z) and a 
temperature field T = Tq + 9. The displacement field 
Ui and the relative temperature field 9, as well as the 
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strain and stress tensors and all functions of 

position and time. 

We consider pure transverse motion Y{x) in the y di- 
rection and make the usual Eulcr-BcrnouUi assumption 
that the transverse dimensions of the beam, b and c, are 
sufficiently small compared with the length L of the beam 
and the radius of curvature R of the bending that any 
plane cross section, initially perpendicular to the axis of 
the beam, remains plane and perpendicular to the neu- 
tral surface during bending. The neutral surface is the 
one running through the length of the beam which suffers 
no extension or contraction during its bending. 

We take the surfaces of the beam to be stress free which 
implies that all but the (Txx component of the stress ten- 
sor vanish on the surface. Because the beam is thin this 
approximately holds in its interior as well. Hooke's law 
for the thermoelastic beam then takes a rather simple 
form 
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Uzx = 0, 
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(9a) 

(9b) 
(9c) 



where E is Young's modulus, a is Poisson's ratio, a is 
the linear thermal expansion coefficient (0), and we have 
taken into account the fact that strain arises both from 
mechanical stress as well as thermal expansion. It is siia- 
ple to show (see, for example. Landau and LifshitzEll) 
that the longitudinal strain component Uxx, a distance y 
away from the neutral surface, is equal to y/R. By re- 
placing the curvature of the beam 1/R with —d'^Y/dx^ 
we may express the non-zero components of the strain 
field in the beam as 
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(10a) 
(10b) 



Following the standard derivation procedure^ for 
isothermal beams with no thermoelastic couplingtiJ but 
with the modified thermoelastic strain ( p^ ) leads to an 
equation of motion for the beam of the form 
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(11) 



where p is the density of the beam, and ^ = 6c is the area 
of its cross section. The quantities / and It are integrals 
over the cross section of the beam giving the mechanical 
and the thermal contributions to its moment of inertia 



I = y dydz 

I A 
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and /t = / yOdydz. (12) 
'a 



In evaluating the moment of inertia I we have neglected 
the deviation of the cross section from its rectangular 
shape which arises from having a non-zero Poisson ratio 



0.5 



0.4 



0.3 



0.2 



0.1 



0.0 







/ \ 

/ 1 (!R(a;)-a;oj 




i ~ 2.225 


/ t^'O 




\ 
/ / 


\ 


1 ^(lo 


/ / 
/ / 
/ / 


\ 

\ 

\ 




/ / 


\ 


\ 


/ / 


\ 




/ / 


\ 






>. 




1 ' 
1 


1 ' 

10 



FIG. 1. Universal plots of the frequency shift and attenua- 
tion, Eq.(p7|), of small flexural vibrations in thin rectangular 
beams due to thermoelastic coupling. 



CT. Such an approximation is justified for small deflec- 
tions since the error it introduces is only on the order of 
the transverse beam dimension divided by the radius of 
curvature of the bending. 

To the equation of motion (|ll ) we add the heat equa- 
tion, which in the presence of thermoelastic coupling is 
given byEll 
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(13) 



We make two simplifications to this equation. First, since 
6 <^Tq, we can safely replace T by Tq in the second term 
on the right-hand side of the equation. Not doing so 
will introduce unnecessary nonlinearities into the prob- 
lem. Second, noting that thermal gradients in the plane 
of the cross section along the j/-direction are much larger 
than gradients along the beam axis, and that no gradi- 
ents exist in the z-direction, we replace by d^9/dy\ 
Substituting the value of the strain field from Eq. ( [l(i| ) 
we finally get 



l + 2Ai 



1 



1-2(7 



89 



dt ^ 9y2 



^Edd^Y 



where we have identified the relaxation strength of 
Young's modulus (^). 



V. SOLUTION OF THE THERMOELASTIC 
EQUATIONS FOR HARMONIC VIBRATIONS 

To calculate the effect of thermoelastic coupling on the 
vibrations of a thin beam we solve the coupled thermoe- 
lastic equations (|l^) and (|l^) for the case of harmonic 
vibrations. We set 
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FIG. 2. Universal plot of thermoelastic damping of small flexural vibrations in thin beams (|29|). The damping is ploted in 
units of the relaxation strength Ae ~ Ea^To/C, as a function of the dimensionless variable ^ = b-^u}o/2x for one decade 
above and below its peak value ~ 0.494, occurring at ^ ~ 2.225. The two bounding Lorentzians (A.3) in the variable ^^/\/24 
are shown along with Zener's approximation Q^^/Ae (^. The relative error in Zener's approximation {Q^^ — Q~^)/Q~^ is 
shown in the inset. 
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calculate the temperature profile along the beam's cross 
section using the heat equation (|l^, and substitute it 
into the equation of motion ( pi] ) in order to obtain the 
normal modes of vibration and their corresponding fre- 
quencies. We expect to find that in general the frequen- 
cies are complex, the real part 3?(a;) giving the new eigen- 
frequencies of the beam in the presence of thermoelastic 
coupling, and the imaginary part |3(w)| giving the at- 
tenuation of the vibration. The amount of thermoelastic 
damping, expressed in terms of the inverse of the quality 
factor, will then be given by 



3(c.) 



(16) 



which is the fraction of energy lost per radian, the factor 
of 2 arising from the fact that the mechanical energy of 
the beam is proportional to the square of its amplitude. 

Substituting ^) into the heat equation ([ij) and ne- 
glecting the term of order on its left-hand side (which 
will only introduce a correction of order A|n to the final 
result) yields the following equation for Oq 
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a dx"^ 



(17) 



whose solution is 
a dx^ 



y — A sm{ky) + B cos{ky) , 



where 



k 




(18) 



(19) 



The coefficients A and B are determined by taking as 
boundary conditions the requirement that there be no 
flow of heat across the boundaries of the beam so that 
d9o/dy = at y = ±6/2. The temperature profile across 
the beam is then given by 
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(20) 



Now that we have the temperature profile we can sub- 
stitute it into the integral It ( |12|) for the cross section's 
thermal moment. Because /, and = Eo^Tq/C 
are all constant along the beam, the beam equation may 
subsequently be expressed as 
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where the complex function /(oj) is given by 
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The equation of motion for the beam (|2|) is formaUy 
identical to that of the isothermal beam with no ther- 
moelastic coupling. The only difference being that the 
isothermal value of Young's modulus E is replaced by a 
frequency dependent modulus 



E,„ = E 



1 



1 + /H 



(23) 



When Lo becomes very large, fioj) 0, and Young's 
modulus tends to its adiabatic, or unrelaxed, value Ead = 
E{1 + Ae)- When lu is very small, /(w) —1, and 
Young's modulus recovers its isothermal value E, as ex- 
pected. For intermediate frequencies E^^ is complex. 

The normal modes of vibration of the beam are given, 
as in the isothermal case, by 

Yo{x) = Asm{qx) + B cos{qx) + C smh{qx) + Dcosh{qx), 

(24) 

where the coefficients A through D and the allowed val- 
ues of q are determined, as usual, by the boundary 
conditions at the two ends of the beam. For exam- 
ple, for beams clamped at both ends or free at both 
ends qnL = a„ = {4.730,7.853,10.996,...}, and for 
cantilevers clamped at one end and free at the other 
qnL = Qn — {1.875,4.694,7.855,...}, in all three cases 
tending for large n to odd-integer multiples of 7r/2. The 
dispersion relation between uj and qn for the thermoelas- 
tic beam is given by 
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(25) 



where wo is the isothermal value of the eigenfrequency. 

Neglecting corrections of order A|; we may replace 
/(w) in the square root by /(wq). The dispersion relation 
(p5|) then becomes 



LO — UJi3 



(26) 



from which we can easily extract the real and the imagi- 
nary parts, giving the thermoelastic corrections of order 
to the eigenfrequencies of the beam as well as the 
corresponding attenuation coefficients 
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(28) 
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FIG. 3. Peak value of thermoelastic damping (5„._ 
~ 0.494A_B (maximum value of the plot in Fig. ^ plotted 
for Gallium Arsenide and for Silicon asrfti fanction of temper- 
ature. Experimentally reported valuefj'E3 are used for the 
thermodynamic quantities E, a, and C. 



The universal behavior of the normalized frequency shift 
(9?(w) — u;o)/ujqAe and of the normalized attenuation 
^{u!)/luoAe as functions of the dimensionless variable ^ 
are shown in Figure 0. 

Using the definition ( p^ of the quality factor we ar- 
rive at an expression for thermoelastic damping in a thin 
beam which, to first order in A^;, is given by 



Ea^Ta 
C 



6 



6 sinh ^ + sin ^ 
^•^ cosh ^ + cos ^ 



(29) 



VI. 



DISCUSSION OF THE RESULTS AND THEIR 
EXPERIMENTAL IMPLICATIONS 



We show in the Appendix that the exact expression 
( p9| ) for thermoelastic damping is tightly bounded be- 
tween two Lorentzians in the variable ^^/v^j and that 
it behaves as for small ^ and as 1/^^ for large val- 
ues of We also compare the exact result with Zener's 
approximate expression (H). The damping is peaked at 
^0 — 2.225 with a maximum value of Q^J^^/A^ ~ 0.494. 
The universal behavior of Q^^/Ae as a function ^ is 
shown in Figure |^. 

The first conclusion to be drawn from this universal 
behavior is that the peak value of thermoelastic damp- 
ing, given approximately by 0.494A^, is independent of 
the dimensions of the beam. It only depends on temper- 
ature through the thermodynamic properties E, a, and 
C of the material. The values of (5~ax ^'^^ GaAs and Sil- 
icon, typical materials used in the fabrication of MEMS 
and NEMS, are plotted as a function of temperature in 
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TABLE II. Thermal relaxation rates i/t = 1/2ttti = /2ttx 
for longitudinal thermoelastic waves in GaAs and Silicon at 



TABLE I. Relaxation strengths Ae{T) and thermal diffu- 
sion lengths £t (in micrometers) for GaAs and Silicon at three 
representative tempcpatures. The values are calculated from 
experimental data^'ES and is used for generating the plots of 
Figure ^ 



Figure ^. Thus, for example, at room temperature one 
expects to observe quality factors no higher than lO'^-lO^ 
if one is operating at around ^ — 2.225. 

The dimensionless variable is proportional to the 
product LOT used in Zener's model of the anelastic solid. 
Because we are mainly concerned vi^ith the dependence 
of thermoelastic damping on scale, and since both w, 
the beam's eigenfrequency, and r, its thermal relaxation 
time, depend on the dimensions of the beam, we find 
it more suitable to explicitly express the dependence of 
^ on the dimensions of the beam. Instead of the usual 
plots of Debye peaks as a function of frequency, we plot 
the damping curve as a function of the dimensions of the 
beam. 

To do that we express the (isothermal) eigenfrequen- 
cies (ESI) in terms of the beam dimensions 




We then get 
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where £t is a thermal diffusion length, proportional to 
the phonon mean free path, given by 



(32) 



We use experimentally reported values i©f the thermo- 
dynamic properties of GaAscJ and Silicone^ to obtain £t 
as well as Ae{T), listed in Table |, for three represen- 
tative temperatures: lOK, lOOK, and 300-?^. We then 
use these values for illustrative purposes to plot the de- 
pendence of thermoelastic damping on geometry in three 
different ways: 

1. Q"^ vs. beam width b for fixed aspect ratio L/b, 

2. Q"^ vs. beam width b for fixed beam length L, 

3. Q"^ vs. beam length L for fixed beam width b. 



values are calculated 
and listed here to il- 



three representative temperatures, 
from experimentally reported data 
lustrate the typical frequencies at which thermoelastic damp- 
ing of longitudinal waves is most significant. The relaxation 
strengths at these temperatures are the same as those listed 
in Table ffl 



The outcome is shown in Figure ^ for the case of a beam 
clamped at both ends vibrating in its fundamental mode 
(a„ = 4.73). The slopes of the curves, plotted in log-log 
scale, are ±1, ±3, and ±2 respectively, which is easily 
understood through the relation (|3^) and the fact that 
grows as for small ^ and decays as ^"^ for large 

t 

It is clear from these examples that thermoelastic 
damping is a significant source of dissipation for MEMS 
and NEMS at temperatures around lOOK and above. 
The reason for this is very simple and follows from the 
fact that as the beam becomes smaller its eigenfrequency 
increases at the same time that its thermal relaxation 
time decreases. The product of the two, which can be 
controlled by independently varying the two dimensions 
b and L, can therefore remain of order unity down to the 
nanometer scale. 



VII. THERMOELASTIC DAMPING OF 
LONGITUDINAL WAVES 

For completeness we would like to say a few words 
regarding thermoelastic damping of longitudinal modes. 
Even though such modes do not come into play when 
considering MEMS and NEMS resonators they may af- 
fect general questions of heat flow and energy relaxation 
in other elements of MEMS and NEMS. The process of 
thermoelastic damping of longitudinal waves is similar to 
that of flexural waves in that energy is dissipated through 
the irreversible flow of heat from hot to cold regions of 
the solid. The difference is that the distance between 
these regions is not fixed by the transverse geometry of 
the device. Because a longitudinal wave is a compression 
wave hot and cold regions are separated by half a wave 
length A along the propagation direction of the wave. 
The wave length A is inversely proportional to the fre- 
quency and therefore the thermal relaxation time r/ for a 
longitudinal thermoelastic wave, which is proportional to 
A^/x, is in fact inversely proportional to the frequency. 
It is therefore the case that high-frequency waves that 
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FIG. 4. Thermoelastic damping in Gallium Arsenide and Silicon thin rectangular beams plotted for diflFerent geometries. In 
all cases the beams are assumed clamped at both ends and vibrating at their fundamental flexural mode. The corresponding 
thermal diffusion lengths £t and relaxation strengths A_B(r), used here, are listed in Table ^ 
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have very short thermal relaxation times are isothermal, 
or relaxed, and low- frequency waves are adiabatic, or un- 
relaxed. This quite counter-intuitive situation is exactly 
the opposite of what one has in the case of flexural vi- 
brations of beams, and is probably the best example of 
a system in which the relaxation time t depends on the 
frequency uj. 

Nevertheless, as can be-seen for example in the work of 
Chadwick and SneddonJlJ thermoelastic damping of lon- 
gitudinal waves still takes the form of a relaxation peak, 
with the isothermal and adiabatic limits exchanged, and 
a with geometry independent characteristic thermal re- 
laxation time Ti — x/cf where q is the speed, or phase 
velocity, of longitudinal waves. Thus, the position of peak 
thermoelastic damping for longitudinal waves is fixed and 
only depends on the thermodynamic properties of the 
material as a function of temperature. Some typical val- 
ues for the thermal relaxation rates of l ong itudinal waves 
in GaAs and Silicon are given in Tabic ||. 



VIII. CONCLUSIONS 

We have established here that thermoelastic damping 
is a significant source of dissipation down to the nanome- 
ter scale. We gave a simple derivation of an exact expres- 
sion ( p9| ) for thermoelastic damping in thin rectangular 
beams, compared this exact expression with Zener's well 
known approximation (j^), and examined the implication 
of our result on micrometer- and nanometer-scale res- 
onators. It is interesting to note, as a consequence of 
our analysis [see Figure^ (a) and (6)], that for beams of 
constant aspect ratio and constant temperature above a 
certain beam width the quality factor increases linearly 
with the size of the beam. This may provide a partial ex- 
planation for the linear increase in dissipation as systems 
become smaller. 

We have made a number of approximations and as- 
sumptions along the way which we would like to summa- 
rize here: 

1. We have derived and solved the thermoelastic equa- 
tions of a thin beam undergoing small flexural vibrations. 
We should not expect our result to strictly hold for beams 
with small aspect ratios L/b. Nevertheless, we do expect 
to see the same kind of behavior, showing a Debye-like 
dissipation peak, even at smaller aspect ratios. We do 
not expect our result to hold for large amplitude vibra- 
tions where the Euler-BernouUi assumption is known to 
fail, and where non-linear behavior begins to take over. 

2. In displaying the expected relaxation strengths 
(Figure 0) and damping curves (Figure § for GaAs 
and Silicon resonators with various geometries we used 
experimentally-reported bulk elastic and thermodynamic 
properties. According to recent moleculer dynamics-|Sim- 
ulations of thin quartz beams by Broughton et alEB one 
is justified in using bulk properties down to beam widths 
of about 10 micrometers. For smaller widths, quantities 



such as Young's modulus change drastically. In any case, 
one should consider the plots only as serving for illustra- 
tive purposes, in particular because some of the thermo- 
dynamic properties, especially below about 200K, vary 
from sample to sample. 

3. The theory of thermoelasticity is valid in a regime 
where thermal phonons arc diffusive and a temperature 
field can be defined locally. One should expect to see 
deviations from this theory when the phonon mean free 
path becomes comparable to the system size, or when the 
relaxation rate of the phonons to their equilibrium Bose 
distribution becomes comparable to the resonator's fre- 
quency. In the first case, where the transport of thermal 
energy crosses over from being diffusive to being ballistic 
one may expect to see thermal relaxation times that are 
linear in the beam width (proportional to b/v where v 
is the phonon velocity) instead of the diffusive quadratic 
dependence 6"/x. The second effect introduces an ad- 
ditional dissipation mechanism which may be viewed as 
caused by the "viscosity" of the phonon gas as it relaxes 
to its equilibrium state. This latter mechanism is called 
the "Akhiezer effect. "EZl We intend to explore these cor- 
rections to the diffusive regime, as well as the fully ballis- 
tic phonon regime, in the context of mesoscopic systems 
in a future pubhcation. 
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APPENDIX A: WHERE IS THE LORENTZIAN? 

Because Zener's approximation ^ using a single 
Lorentzian in the variable ut^ is so good, it is illumi- 
nating to tryto identify this Lorentzian behavior in our 
own result (29). To do so let us expand both numerator 
and denominator oi / Ae in powers of ^, 



Q- 



^ 6 



.■SI'S 



Ac2 



2-5! 
9! 



13! 



1 



12! 



(Al) 



One can clearly see that for small values of ^, the two 
series in parentheses each tend to 1, their ratio approach- 
ing 1 from above as ^ decreases. For large values of ^, 
the leading 1 in the denominator may be neglected, and 
one can show with a little bit of effort that the ratio of 
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the two series in parentheses tends to 5/4 from below as 
^ increases. Defining £ as the Lorentzian 



V 



(A2) 



we see that for any value oi Q ^/Ae is bounded be- 
tween the two Lorentzians 



< 01 < ^c(J^ 

- Ae - 2 V\/24 



(A3) 



Figure |^ shows the universal thermoelastic damp- 
ing curve Q~^/Ae in relation to these two bounding 
Lorentzians. The inset shows the difference between 
Zener's Lorentzian approximation, which in the above 
notation takes the form 



(A4) 



and the exact result (|9|). Note that ~ 4.935 and 

^/24 ~ 4.899 differ by less than 1%. It should be em- 
phasized that on the isothermal side of the peak (low 
frequencies) the two expressions differ by less than 2%, 
which is the error anticipated by Zener in keeping only 
the first term in his expansion. On the adiabatic side of 
the peak (high frequencies) the error increases in the first 
decade to as much as 15%, reaching 20% in the limit of 
infinite ^. 
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